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Abstract. A possible mechanism for the formation of a 90° misfit dislocation at the 
Ge/Si(001) interface through homogeneous nucleation is identified from atomic scale 
calculations where a minimum energy path connecting the coherent epitaxial state 
and a final state with a 90° misfit dislocation is found using the nudged elastic band 
method. The initial path is generated using a repulsive bias activation procedure 
in a model system including 75000 atoms. The energy along the path exhibits two 
maxima in the energy. The first maximum occurs as a 60° dislocation nucleates. The 
intermediate minimum corresponds to an extended 60° dislocation. The subsequent 
energy maximum occurs as a second 60° dislocation nucleates in a complementary, 
mirror glide plane, simultaneously starting from the surface and from the first 60° 
dislocation. The activation energy of the nucleation of the second dislocation is 
30% lower than that of the first one showing that the formation of the second 60° 
dislocation is aided by the presence of the first one. The simulations represent a step 
towards unraveling the formation mechanism of 90° dislocations, an important issue in 
the design of growth procedures for strain released Ge overlayers on Si(100) surfaces, 
and more generally illustrate an approach that can be used to gain insight into the 
mechanism of complex nucleation paths of extended defects in solids. 
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1. Introduction 

Heteroepitaxial systems play an important role in modern microelectronics technology 
[lj. Due to lattice mismatch between the film and the substrate, considerable elastic 
strain energy is accumulated during epitaxial growth. For sufficiently thick films, 
the strain energy is released through the formation of defects, leading to the loss of 
coherent epitaxy Controlling the film-substrate interface quality and preventing 

defect formation within the film is an important goal in modern technology. This 
requires detailed information about the complex atomic rearrangements during strain 
relaxation. Current experimental methods do not, however, allow the measument of 
detailed atomistic evolution in such processes. Therefore, theoretical modeling can play 
an important role in studies of defect formation mechanisms and atomic structure. 

Misfit dislocations (MDs) are the most important type of defects in relaxed 
heteroepitaxial systems [3]. Early theoretical studies of MDs were based on the 
comparison of the energy of configurations with and without MDs P2E!. Later on, it has 
been recognized that the formation of defects leading to the loss of coherent epitaxy is a 
thermally activated process, as illustrated by the experimentally observed temperature 
dependence of the critical thickness of the film [2], The state of the epitaxial system 
formed in experiments is thus not only determined by the nature and energy of the 
final state, but also by the kinetic factors which depend on the transition path, and in 
particular on the energy barriers that need to be overcome to achieve strain relaxation. 

The Ge/Si heterostructure is a particularly important heteroepitaxial system in 
microelectronics applications. It is used in optoelectronics and for creating high speed 
electronic components mom. It is also an important substrate for GaAs growth 
on silicon pm ITT] . The Ge lattice constant is around 4% larger than that of Si, 
so a perfect, coherent film can only grow up to a few monolayers. There are two 
competing mechanisms for relaxation of the him strain in Ge/Si(001), namely through 
MD nucleation, and by 3D island formation in the Stranski-Krastanow growth mode. 
Experiments show that the Stranski-Krastanow growth mode can be suppressed by 
proper choice of growth conditions (for example by using low temperature [12] and by 
using surfactants 0)- There are two main types of MDs in the diamond structure, 
namely 60° and 90° MDs. The terminology reflects the angle between the direction of 
the Burgers vector and the dislocation line. The nucleation of 60° MDs in pure Si near a 
surface step under external stress has been extensively studied [13] dU d5] H6] HTJ dHl [19] . 
It has been shown that since the 60° MDs can glide on the dense (111) planes, they can 
form through half-loop nucleation on these planes. 

The Burgers vector of a 90° MD lies in the (001) plane while the Burgers vector of 
a 60° MD does not. A 90° MD thus releases larger misfit strain than a 60° MD. The 90° 
MD, however, is sessile, i.e. it is relatively immobile. Experimentally, a regular array of 
edge MDs lying at the interface and almost entirely releasing the strain from the film 
has been obtained [201 El] • Since 90° MDs are sessile, the 90° MDs must form through 
the reaction of non-sessile dislocations. Bolkhovityanov et al. have presented a review 
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of various possible mechanisms for the formation of 90° MDs [21] • All the mechanisms 
described involve the reaction of two complementary 60° MDs. Such a reaction could 
for example be of the form: 

a/2[011](lll) + a/2[10l](lll) ->• a/2[110](001), (1) 

where a is the Si lattice constant, a/2 [Oil] (111) indicates that the MD has a Burgers 
vector of a/2[011] and that it glides on the (111) plane. In the induced nucleation 
mechanism, the presence of a 60° MD favors the formation of a complementary 60° 
MD on the mirror plane and the two MDs react. Two ways for the complementary 
MD to nucleate have been considered. The complementary MD can either form by a 
half-loop nucleation from the surface [22] [23., 24] or it can form from the existing MD 
[25] following the reaction: 

a/2[011](lll) ->• a/2[T01](lll) + a/2[110](001). (2) 

To the best of our knowledge, only one atomistic simulation of the nucleation of a 
90° MD through the induced nucleation of a complementary 60° MD has been presented 
so far pEj. However, a simplified quasi two-dimensional model of the Ge/Si(100) 
system was used and surface reconstruction was not included, so the activation energy 
for nucleation, for example, is not expected to be accurately estimated. To better 
understand the process of the 90° MD formation, one has to find minimum energy paths 
for this transition on the energy surface in a large enough atomic scale representation 
of the system. This is a nontrivial problem due to the high dimensionality of the 
configuration space. 

We have recently explored the atomic relaxation mechanisms at the microscopic 
scale in three-dimensional heteroepitaxial systems with hexagonal symmetry [27]. 
Starting from the initial epitaxial state, we generated the final configuration containing 
various kinds of defects with the repulsive bias potential method (RBP) [28]. Then, 
the nudged elastic band method [29l ; 30] was employed to determine a minimum energy 
path for the transition. This approach allowed us to classify the atomic mechanisms 
leading to strain relaxation in two and three dimensional systems. 

In the present work we apply this procedure to a more complex system, an atomistic 
model of the Ge/Si(100) system including the (2 x 1) surface reconstruction. We 
note that under experimental conditions where defects and impurities are present, the 
nucleation of dislocations may not occur through the homogenous nucleation process. 
Due to lack of detailed experimental data realistic modeling of such processes is often 
challenging. In the present work we consider homogenous nucleation only for the sake 
of conceptual clarity. A minimum energy path for the formation of a 90° MD is found. 
A possible mechanism is identified that starts with the nucleation of a 60° MD. Then, 
a complementary 60° MD nucleates with a lower activation energy and the two MDs 
finally react to form the 90° MD. While the activation energy is estimated to be very 
high, too high to be overcome by thermal activation alone, the mechanism identified 
illustrates how the presence of one MD can induce the formation of a second one. The 
origin of the high activation energy is probably related to limitations in the model, such 
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as the form of the interaction potential used and the absence of defects in the initial, 
epitaxial state. Further work is needed to resolve these issues. 

2. Model 

The computational model system has the shape of a parallelepiped, as shown in Fig. 
ED with periodic boundary conditions applied in the x — y plane to mimic an extended 
slab consisting of 31 atomic layers of Si (representing the substrate), and 19 layers of 
Ge (representing the him). The crystal orientation of the surface of the substrate is 
(001). The x, y and z axes are oriented along the [110] , [110], and [001] directions, 
respectively. The bottom two layers of the substrate were fixed to mimic interaction 
with an extended Si crystal below the slab and to prevent the system from moving as a 
whole. 

Since dislocations create long distance strain fields, it is important to use as large 
a system as possible to minimize finite-size effects. We checked that the system is large 
enough in the sense that further increase in the substrate thickness or in the system 
size in the x and the y directions does not significantly change the energetics of the 
transition paths calculated here. Here, the aim is to generate dislocation lines oriented 
along the y axis ([110] direction). Because of periodic boundary conditions, a single MD 
in the computational cell corresponds to a network of parallel MDs. From experimental 
observations, the average distance between such MDs has been measured to be around 
100 A [20]. Therefore, we chose the size of the system in the x direction to accomodate 
50 atomic rows. The lateral dimension of the system (in the x — y plane) is thus 50 x 30 
atomic rows, corresponding to a system of size of 19.2 x 11. 5 x 6.8 urn 3 . The total 
number of atoms in the model is 75000. 

The canonical Stillinger-Weber (SW) potential function is used to describe the 
interatomic forces [31]. The SW potential is, of course, only a crude approximation 
to the atomic interactions, but it is the most widely used potential function and we 
have therefore adopted it here. The system size needed to describe the phenomena 
studied here makes the computational effort in quantum mechanical density functional 
calculations prohibitively large. The SW potential enforces sp 3 hybridization while 
in reality the bonding arrangements are more flexible [32] . As a consequence, the 
results presented here are expected to give an overestimate of the activation energy 
for the dislocation formation. The (001) surface is dimer reconstructed (cf. Fig. 1) 
in the calculations, as has been observed during growth of Ge films on Si(001) [ 33] . 
Experimental observations, furthermore, show a complex 2 x N reconstruction pattern 
on the Ge/Si(001) surface where every N -th dimer is missing. We used the ideal (2 x 1) 
reconstruction as an initial state in our model for the sake of simplicity. The parameters 
of the Stillinger-Weber potential are given in Table [0 They are the same as the one 
used in [Ml SI] • 


5 


Table 1 . Parameters of the Stillinger-Weber potential for mixed Ge-Si 
interactions,taken from ref. |34j . 


Parameter 

Si-Si 

Ge-Si 

Ge-Ge 

cr 

2.0916431 

2.1354 

2.17912051 

e 

2.17 

2.0427 

1.93 

A 

21.0 

- 

31.0 


3. Method 

Epitaxial films with thickness beyond the equilibrium critical thickness will tend to 
relax to a state that includes defects, but this process typically inolves overcoming an 
energy barrer A E during the transition from the coherent state. To find the activation 
energy and determine the mechanism of the transition, a minimum energy path for 
the transition needs to be found. A direct simulation of the classical dynamics of the 
atoms is typically not useful because the transition is a rare event on the time scale 
of the atomic vibrations. While the transition rate can be increased by increasing the 
temperature, a cross-over to a different mechanism, favored by entropy, is likely to occur. 
To clarify the mechanism of the relaxation, the identification of a minimum energy path 
is more useful, ffowever, this is a challenging task as many atoms are involved and it is 
not clear what the atomic structure of the final, relaxed state is. 

To generate an isolated defect at a given location in the sample, we have previously 
introduced [ 28 j a particularly simple but efficient method called the repulsive bias 
potential (RBP) method. In the RBP method, the system is placed in a fixed, external 
repulsive bias potential which makes the initial state unstable 

U to t{r,f 0 ) = U(f) + Aexp{-[(f- f 0 )/a] 2 }. (3) 

ffere the components of r 0 and r contain all the atom coordinates of the initial and 
current configurations, respectively. The U(r) term contains the atomic interactions 
(here the Stillinger-Weber potential function) and the RPB potential is the additional, 
spherically symmetric Gaussian term with strength A, range a and a maximum value 
at r 0 . When A and a have been chosen appropriately, forces computed from Eq. (J3)) 
can displace the system from its initial state to a nearby local minimum on the energy 
surface. In practice, this is done by applying an energy minimization using f/ to t as 
objective function. 

First, the initial epitaxial state is prepared by minimizing the total energy. We have 
used a procedure based on classical dynamics where the velocity of each atom is set to 
zero whenever it has a component opposite to the direction of the acceleration |29l 130] . 
The standard leap-frog algorithm was used to numerically integrate the equations of 
motion. The minimization was considered to be converged when the maximal force 
acting on an atom had dropped to less than 10 -4 eV/A. 
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To generate an atomic configuration for the final, defected state, the RBP bias is 
turned on and a group of atoms is displaced from the initial position to bring the system 
closer to the final state of the particular relaxation process we want to study. Rather 
than trying random initial displacements, some knowledge of the defect generation 
mechanism is useful for expediting the process. Then, the total energy minimization is 
reapplied to determine the structure at the new, defected local minimum. 

It is important to note that this method can generate many different final states 
depending on both the choice of initial displacements and the exact form of the RBP. 
By making the repulsive bias sufficiently localized around the initial potential minimum, 
the final state energy depends only on the true potential of the system and not on the 
fictitious repulsive bias. We consider only final configurations with precisely one isolated 
defect, namely a 90° MD to unravel the mechanism and energetics of its formation. 

To generate an isolated 90° MD at certain position in the system, one needs to use 
preliminary knowledge about the geometry of the diamond lattice and the arrangement 
of atoms in the dislocation core to select particular set of displacements and parameters 
of the RBP. Our trial and error attempts resulted in the following procedure: Atoms in a 
wedge-shaped triangular prismoid were displaced by 1.5 A in the direction of [Til]. The 
RBP parameters were chosen to be A = 2000 eV and a = 0.011/A 2 . The minimization 
of [/ to t then suffices to move the system away from the potential basin corresponding to 
the perfect, epitaxial state to the energy minimum corresponding to a 90° MD state. 

With both an initial and final state atomic configuration available, the nudged 
elastic band (NEB) method [29j, [30], is used to find a minimum energy path (MEP) 
between the two. A set of intermediate configurations (’images’ of the system) are first 
generated to create a discretized path between the initial and final states. Iterative 
displacements of the images along the NEB force using some minimization algorithm 
then gradually bring the images to an MEP. When more than one MEP on the energy 
surface connects the given endpoints, the NEB minimization typically converges to the 
MEP closest to the initial path. The initial location of the images is often obtained by 
linear interpolation of the atomic coordinates between the initial and final states. For 
the present calculations however, we found that this can lead to numerical instabilities 
due to the strong hard core repulsion of the interatomic potentials as atoms are brought 
close to each other. A method involving interpolation of pairwise distances has recently 
been presented to circumvent this problem [35]. Here, however, we used the intermediate 
configurations obtained during the minimization with the repulsive bias. This approach 
leads to fast convergence of the NEB without the instabilities encountered with the 
linear interpolation scheme. The same procedure has previously been used to study 
MD nucleation in metal-on-metal systems with the FCC(lll) surface orientation [2?j . 
In the NEB optimization, we use the same procedure based on classical dynamics with 
velocity zeroing as for the RBP minimization [29, 130]. To get reasonable resolution of 
the transition path, we used about 130 images in the NEB. 


7 


4. Results 

We will first discuss the influence of the him thickness on the energy of the states with 
either a 60° MD or a 90° MD as compared to the energy of the coherent epitaxial state. 
Then the MEP found for the formation of the 90° MD is described. 

4-1- Relaxation energy 

Starting from the coherent epitaxial state (Fig. [2b), the formation of a 60° MD (Fig. 
Eb) or a 90° MD (Fig. [2b) decreases the him strain by moving atoms out of the surface 
layer to form an island of adatoms. For a 90° dislocation, the island is hat and its width 
is proportional to the him thickness. The energy difference A E g = — E\ between the 

coherent, epitaxial state and a state where a MD has formed arises from competition 
between decrease in energy due to strain release in the him and increase in energy 
associated with the formation of the dislocation core and surface defects. To study this, 
the Ge him thickness was varied between 9 and 19 ML while keeping 31 atomic layers in 
the Si substrate. Fig. [3] shows how A E g corresponding to the formation of a 60° MD and 
of a 90° MD vary approximately linearly with the thickness of the Ge him. This linear 
behavior arises from the fact that the decrease in strain energy is almost proportional 
to the him thickness while the energy of dislocation core and surface defects energy is 
independent of the him thickness for films thicker than a few MLs. The results show that 
the energy of the system is lowered by forming a 60° MD beyond a critical him thickness 
of 19 atomic layers, which corresponds to approximately 2.7 nm. This estimate agrees 
well with the results of Ichimura et al. [26]. The energy of the system is lowered by 
forming a 90° MD beyond a critical him thickness of 10 atomic layers, which corresponds 
approximately to 1.5 nm. This estimate agrees well with experimental observations [36], 
but is a little higher than that obtained by of Ichimura et al, 0.8 nm [26] . 

4-2. Minimum energy path 

The energy along the MEP found for the formation of a single, straight 90° MD at 
the interface starting from the coherent epitaxial state is shown in Fig. [U A few 
intermediate configurations along the reaction path are shown in Figs. [5] and [6] A video 
can be downloaded from the supplemental material 0. 

The reaction can be divided into two parts. The hrst part corresponds to the 
formation of a 60° MD with Burgers vector [Oil] and has an activation energy of 54 
eV. The Burgers vector was determined using the dislocation extraction algorithm [37J. 
Figures [5b and [6b represent the configuration of the saddle point for the formation of 
the 60° MD. They clearly show that the 60° MD forms through a half-loop nucleation on 
the (111) slip plane. At the saddle point the dislocation has reached the interface and 
can be described as a straight dislocation at the interface terminated by two threading 
arms reaching the surface. The dislocation then spreads by glide of the two threading 

f http://figshare.com/articles/Movie_ogg/1423292 



arms in opposite directions. This dislocation growth is associated with a decrease of the 
energy due to the strain release in the film. The energy decrease is small in this case 
because the film thickness corresponds to the critical thickness for the formation of a 60° 
MD. Due to the periodic boundary conditions, the threading arms meet (see Figs. [5jo 
andEb) and annihilate each other leaving only a straight 60° MD at the Ge/Si interface 
(shown in Figs. EF andEF)- This annihilation leads to a decrease of the dislocation 
length and a decrease in energy. The straight 60° MD corresponds to the intermediate 
local minimum in the reaction path and has slightly lower energy than the coherent 
state. Figure [2b shows that the formation of the 60° MD also induces the formation of 
a double layer step on the surface. 

In the second part of the reaction, a complementary 60° MD with Burgers vector 
[10T| forms in the (111) plane which is the mirror plane of the (111) [§|. The activation 
energy for the formation of this complementary MD is around 37 eV, 30% lower than 
the activation energy for the formation of the first 60° MD. The presence of the first 
MD thus assists the formation of the complementary MD on the mirror plane. 

In previous studies, it has been assumed that the complementary MD nucleates 
either from the free surface [21] or from the existing MD [22] • In our calculations, both 
events occur simultaneously as one loop forms from the existing MD while a half-loop 
forms from the surface (see Figs. [5bl andEFO- The two loops then merge to leave only 
the threading arms shown in Figs. EF andEF- The threading arms can then glide and 
react with the straight 60° MD to form a 90° MD according to Eq. [T] 

The final state of the system with a single 90° dislocation is shown in Figs. EF, EF 
and |2b It is characterized by the appearance of an extended double layer island on the 
surface running along the y direction. 

5. Discussion 

The estimate of the activation energy obtained here, 54 eV, is obviously very high. An 
event with such a high activation energy would never occur by thermal fluctuations 
even if the temperature were close to the melting temperature. However, it has been 
experimentally shown that the half-loop nucleation of a 60° dislocation can occur even 
at a lower strain in a Ge 0 .32Sio.68/Si(001) film[40]. The mechanism presented here 
might therefore be possible. This implies that we are significantly overestimating the 
activation energy in our calculations. Unfortunately no reliable experimental estimate 
of the activation energy is available in the literature. 

Many atomistic calculations of dislocation nucleation in semi-conductor materials 
have predicted a very high activation energy [n m sn E]. Below, we discuss possible 
explanations for these large values. 

§ According to |24j . the complementary MD is more likely to form not on the exact mirror plane but 
on a slightly translated mirror plane. After a glide of the complementary MD, a pair of coupled 60° 
MDs can be observed near the interface [35] [32] • In our calculations, a 90° MD is present in the final 
state and the complementary MD thus has to form on the mirror plane. 
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5.1. Accuracy of the potential 

The SW potential enforces chemical bonding according to sp 3 hybridization while in 
reality the bonding arrangements can be more flexible [32] . The chemical bonding 
between atoms in the core of the dislocation is, in particular, expected to deviate from 
sp 3 hybridization and thereby be overestimated in the present simulations. In order to 
obtain an estimate of this effect, the difference in energy per atom was calculated for 
each atom to identify which region of the dislocation is primarily responsible for the 
activation energy of the formation of a 60° MD. The total increase in energy of the core 
atoms was found to be 61.7 eV. If the energy of core atoms is overestimated by 30%, the 
activation energy would be overestimated by around 20 eV. In principle, our calculations 
could be corrected by estimating the error in the dislocation core energy from comparison 
with more accurate density function theory (DFT) calculations. However, estimating 
dislocation core energy from DFT is a challenging task [42] 43] [44] and to the best of 
our knowledge no estimate of dislocation core energy in Ge is available. 

Preliminary calculations have also been carried out using a different potential 
function, the Tersoff potential ga m szi, and the activation energy is then found 
to be 12 eV lower than with the SW potential. The large difference in the activation 
energy calculated with these two potential functions indicates that a more accurate 
description of the atomic interactions allowing for deviations from sp 3 hybridization 
could give a significantly lower activation energy. Other potentials such as the Lenosky 
potential [48] have been shown to more accurately decribe dislocation cores in Si [09] 
but a parametrization for Ge is not available. There is clearly a need for developing a 
more accurate potential function for modeling dislocation cores in the Ge/Si system. 

5.2. Absence of defects 

Dislocations usually nucleate from pre-existing defects. For instance, it was found 
experimentally that most of the 60° MD half-loops nucleated from the same location 
in a Geo.32Sio.68/Si(001) him [40]. This observation was interpreted as indicating that 
a defect was initiating the formation of the dislocation. Such defects can be impurities 
[501 HU US] or, more importantly, steps on the surface. When starting from a defect free 
him, the formation of MDs leads to the formation of double-layer steps on the surface 
as shown in Fig. [2] This is energetically unfavorable since step atoms have dangling 
bonds. At the saddle point the increase in energy due to the step edge atoms is 7.7 
eV. If the him initially contains steps as indicated in Fig. [7(a), the formation of an 
MD eliminates the surface steps, as illustrated in Fig. [7(b). The activation energy for 
forming an MD would, thereby, be lowered by the presence of steps in the initial state. 
Steps, furthermore, act as stress concentrators and most of the atomic scale simulations 
of the formation of 60° MD in diamond systems have, indeed, been carried out for 
systems that include a surface step in the initial state HZl HH1H91H31HH H61 EH [52]. 

Since at the saddle point, the increase in energy due to the step edge atoms is 7.7 
eV, we can expect that having a straight step in the initial configuration would lower 
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the activation energy by about 15 eV. The lowering of the activation energy of the step 
should be similar for the nucleation of the initial and of the complementary MD but 
a significant difference should be noted. Formation of an initial MD will be favored in 
the vicinity of a step whereas in order to favor the formation of a complementary MD a 
step has to be located on the mirror (111) plane of an existing 60° MD. During growth 
of the film, it can be expected that steps are propagating along the surface as atoms get 
deposited on the surface and that a step will at some time be located at just the right 
position for inducing the nucleation of the complementary MD. 

6. Summary and Conclusions 

Misfit strain relaxation in the heteroepitaxial Ge/Si(001) system has been studied using 
an atomistic simulation method where a repulsive bias potential is used to generate 
atomic coordinates of a 90° MD and a nudged elastic band calculation is then carried 
to to find a minimum energy path between this and the coherent epitaxial state. The 
Stillinger-Weber interatomic potential energy function has been used to approximate the 
atomic interactions. The minimum energy path found here for the formation of the 90° 
MD involves half-loop nucleation of a 60° MD and the subsequent induced nucleation 
of a complementary 60° MD. The activation energy for the formation of the first and 
second 60° MD is 54 eV and 37 eV, respectively. The presence of the first MD thus 
assists the formation of the second by lowering the activation energy by 30%. The 
large value of the activation energy obtained here may be attributed to the absence of 
defects in the initial state of the him and/or to the inflexibility of the Stillinger-Weber 
potential function which does not account for bonding arrangements that deviate from 
sp 3 hybridization and thus gives very high energy for the atoms in the dislocation core 
of a 60° MD. 
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z=[001] 


x=[110] 


Figure 1. The model of the Ge/Si(001) system in the inital, coherent state. 
Blue spheres represent the Si atoms and red represent the Ge atoms. The surface 
is reconstructed as rows of dimer form between surface atoms. 



Figure 2. (a) Side view of the initial, coherent state, (b) Side view of the straight 

60° misfit dislocation, (c) Side view of the straight 90° misfit dislocation. Blue spheres 
represent the Si atoms and red spheres the Ge atoms. 
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Figure 3. Energy difference due to formation of a 60° (green square) and a 90° (blue 
dot) misfit dislocation (normalized to the length of the dislocation line) as function of 
the Ge film thickness. Dotted lines show a linear fit to the data. 
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Figure 4. Energy along the minimum energy path between a defect free, coherent 
initial state and a final state with a single 90° dislocation at the Ge/Si interface. The 
configurations corresponding to the dashed vertical lines are shown in Figs. [5] and [6] 
The Ge film contains 19 atomic layers while the Si(001) substrate contains 31 atomic 
layers. 
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Figure 5. Six intermediate configurations along the minimum energy path for the 
formation of a 90° dislocation at the Ge/Si interface. The position of the configurations 
in the energy profile is indicated by vertical dashed lines on Fig. [4] (a) corresponds to 
the first saddle point, (b) corresponds to the closing of the half-loop due to periodic 
boundary conditions, (c) corresponds to a straight 60° dislocation, (d) corresponds 
to the second saddle point, (e) illustrates the growth of the 90° dislocation, and (f) 
corresponds the straight 90° dislocation. Only atoms whose surrounding does not 
correspond to a diamond lattice according to a common neighbor analysis [53} 154] as 
implemented in OVITO are shown [55] . The (111) and (111) planes are shown in red 
and blue. The (111) plane is not shown in figures (d) and (e). The location of the 
Ge-Si interface is indicated by the gray plane. 
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Figure 6. Schematic representation of the configurations shown in Fig. [5] illustrating 
the mechanism for the formation of 90° misfit dislocation at the Ge/Si interface. The 
light green, dark green and gray lines indicate 60° dislocations with Burgers vector 
[Oil], [101] and [T01], respectively. The white line indicates a 90° misfit dislocation 
with Burgers vector [110]. The (111) and (111) planes are shown in red and blue. 
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Figure 7. Side view of two configurations for a Ge/Si(001) film. Starting from a 
film containing steps (a), the formation of a straight 90° MD leads to the elimination 
of the steps (b). 






























